function v = symlimit(f, var, a, dir)
    % 求表达式的极限
    % v = symlimit(f, var, a) 求表达式f中的变量var在a点的（左右）极限
    % v = symlimit(f, var, a, dir) 求表达式f中的变量var在a点的方向极限，dir可以是"+","-","+-"。
    if nargin < 3; error("输入参数不够"); end;
    if nargin < 4; dir = "+-"; end
    if nargin < 3; a = 0; end
    if any(strcmp(["left", "-"], dir))
        dir = "-";
    elseif any(strcmp(["right", "+"], dir))
        dir = "+";
    else
        dir = "+-";
    end
    sp = sympy_sp();
    v = sp.limit(f, var, a, dir).simplify();
end